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Abstract 

Structure of Mach cones in a crystalline complex plasma has been studied experimentally using 
Q-t an intensity sensitive imaging, which resolved particle motion in three dimensions. This revealed a 

O \ previously unknown out-of-plane cone structure, which appeared due to excitation of the vertical 

wave mode. The complex plasma consisted of micron sized particles forming a monolayer in 
a plasma sheath of a gas discharge. Fast particles, spontaneously moving under the monolayer, 
created Mach cones with multiple structures. The in-plane cone structure was due to compressional 

> 

Q\ [ and shear lattice waves. 
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Complex (or dusty) plasmas are weakly ionised gases containing dust. Due to absorp- 



tion of ambient e 



ectrons and ions of the plasma, dust particles acquire significant negative 



electric charges [l-^. They interact strongly with the plasma and with each other and they 
can form strongly coupled systems analogous to colloids. Dust particles are often confined 
in the sheath region of electrical discharges, where the electrostatic force is strong enough 
to compensate their weight. The sheath region confines the particles strongly in the ver- 
tical direction and they can form monolayer crystals 

aa. 

These systems extend in the 

'third' vertical direction to their finite width which depends on the strength of vertical con- 
finement pj. A vertical "out-of-plane" particle motion is allowed making the dust system 
geometry three-dimensional. However, the magnitude of the vertical displacement is very 
small compared to the interparticle distance (tens of micrometers compared to hundred of 
micrometers). For this reason, monolayer complex plasma crystals are often referred as two- 
dimensional or quasi-two-dimensional (quasi-2D) complex plasma crystals js]. A monolayer 
plasma crystal has not to be confused with a dust multilayer system 9j. 

If a body moves in a wave sustaining medium, it creates a disturbance, also known as a 
wake. V-shaped wakes, also known as Mach cones, are created behind a body moving super- 
sonically. These cones are observed behind boats, supersonic jets, in fluid filled boreholes, 
and in complex plasmas. The structure of Mach cones is determined by the wave modes 
that exist in the medium, with each mode having a distinct contribution. Quasi-2D lattices 
sustain two in-plane wave modes with acoustic dispersion. One of them is compressional 



10(. Compressional cones 



rave been observed 



(longitudinal), the other is shear (transverse) 
in a quasi-2D complex plasma using spontaneously accelerated particles llMl3|. They often 
have a multiple structure with a c omp ressional cone followed by a rarefactional and then 
possibly by another compressional 14J. The cone angle /i is determined by the Mach cone 
rule sin \i = c/V, where V the speed of the supersonic particle and c is the speed of the 
wave that comprises the cone. If the waves are dispersive (e.g. water waves) or change 
their speed as they propagate in an inhomogeneous medium, the Mach cone angle will also 
change [ljj] . Wakes produced by shear waves have been observed experimentally using laser 
excitation [if], [l7]. Since the shear mode is about 4-5 times slower than the compressional 
one, the shear cone fits within the compressional cones, so that both structures can be 
observed simultaneously. Shear cones have been generated by supersonic dislocations in 
stressed complex plasma crystals 



Since the strength of the vertical confinement in quasi-2D complex plasma systems is 
finite, there is a third fundamental wave mode: vertical transverse wave (VTW). It is asso- 
ciated with propagating out-of-plane or vertical oscillations [s], Q, 20] and it has an optical 



dispersion. It is then natural to ask whether a vertical Mach cone can be produced by VTW. 
Optical wave packets propagate with the group velocity, which is, in case of VTW, usually 
slightly lower than the shear wave speed Thus it should be possible to observe com- 
pressional, shear and vertical cone structures at the same time. Here we report a study of 
the 3D structure of Mach cones using an intensity sensitive imaging technique % ]8| , which 
measures vertical displacements of particles in a monolayer lattice. This method has allowed 
us to resolve compressional, shear, and vertical cone structures. 

The experiments were performed in two very similar capacitively coupled radio frequency 



,8 



at Max Planck Institute for Extraterres- 



(RF) glow discharge chambers at 13.56 MHz, one 
trial Physics in Garching (experiment I) and the other [2j| at the department of Electrical 
Engineering and Electronics at the University of Liverpool (experiment II). The experi- 
mental parameters are listed in Table [H The chambers were initially pumped down using 
turbo-molecular pumps. The working gas pressure was maintained by a small argon flow set 
to minimise the disturbance to the lattice. The RF power was applied between the lower 
electrode (^200 mm in diameter) and the grounded chamber. Due to electrode asymmetry, 
the powered electrode acquired a negative self-bias voltage. Typical plasma parameters have 
been measured with a Langmuir probe in the bulk discharge at a pressure of 0.66 Pa and 



discharge power of 20 W [22| yielding the electron temperature T e = 2.5 eV and the electron 
density n e = 2 x 10 9 cm" 3 . A monolayer particle suspension (50-60 mm in diameter) was 
formed by injecting melamine formaldehyde microspheres (9.19±0.14 nm in diameter, with 
a mass of 6.1xl0~ 13 kg) into the plasma sheath above the lower electrode. 

The microparticles were illuminated by a thin horizontal laser sheet, which had a Gaussian 
profile in the vertical direction and a uniform profile in the horizontal direction. The sheet 
thickness was nearly constant across the crystal. The particles were imaged through the 
top window by a Photron FASTCAM high speed camera. The horizontal coordinates x and 
y as well as the velocity components v x and v y of individual particles were then measured 



with sub-pixel resolution using standard particle tracking techniques 231-4251] . An additional 
side-view camera was used to verify that our experiments were carried out with a single 
main layer of particles. In order to determine the relative vertical positions z and velocities 



TABLE I. Parameters used in the experiments and simulation. 



Parameter 


Symbol 


Unit 


Experiment I Experiment II Simulati 


argon pressure 


P 


(Pa) 


u.00 


1 A^ 
1.410 




RF power 


P 


(W) 




Q 
O 




camera frame rate 




(fps) 


9^n 
zou 


IOC 

1Z0 




average interparticle distance 


a 


(jtm) 


000 


0(0 


AO^ 
4Z0 


longitudinal dust-lattice wave 


speed Cl 


(mm/s) 


01. i 


0"i A 


*}R 1 
oo. 1 


screening parameter 


k = a/Xo 




1 9"? 
1 .ZO 


1 /in 

1 .4U 


n Q/i 


Debye length 


Xd 


(H 


A ^9 


a 1 n 

4tlU 


/I ^9 
40Z 


particle charge in electron charges Z 


00 


loUUU 


i a snn 

14oUU 


1 snnn 

loUUU 


vertical resonance frequency 


fv 


(Hz) 


26.0 


16.5 


23.9 


speed of perturbing particle 


V 


(mm/s) 


40.2 


19.2 


40.0 


longitudinal Mach number 


M L = V/C L 




1.2 


0.8 


1.05 


longitudinal Mach angle 


Hl = arcsin(l/M^) 


(deg.) 


56 




72.5 


measured horizontal damping 


lh 


(s- 1 ) 


0.83 ±0.21 


1.57 ±0.63 




theoretical damping (Epstein) 


1th 


(s- 1 ) 


0.79 


1.73 


1.0 



v z of individual particles, the laser sheet was set slightly above the particle layer, so that the 
position of the intensity maximum was 100 /mi to 200 /im higher than the average levitation 
height of the grains. The intensity of light scattered by the particles depended on their 
vertical displacement from the average height. The vertical velocity was calculated from the 
change of particle intensity in two consecutive frames 7, s|. 

The complex plasma lattices used in our experiments have been characterised by mea- 
suring their main parameters, which are shown in Table [I] The charge and the screening 
parameter in both experiments were determined from the long wavelength phonon spectra of 
the thermally excited lattice oscillations. The spectra have been analysed yielding the lon- 



gitudinal and transverse dust-lattice wave speeds with the method of Ref . 26( . The particle 



charge and the screening parameter were then calculated using the formulae of Ref. [27| . The 
interparticle distances were measured directly from the lattice images, and the Debye length 
was calculated. We have also measured the vertical resonance frequency, which characterises 
the strength of the vertical confinement. 



Mach cones were created using the experimental technique of Refs. Ill, [l2| . A few heavy 
particles (which are often found in many experiments) moved under the main lattice layer. 
They were spontaneously accelerated, most likely by the wake fields of particles in the main 



layer 



281 ]. if the RF power was high enough. Their speeds were nearly constant, limited by 



the neutral gas friction (Epstein drag), and their trajectories were most often straight lines. 
They produced Mach cones due to electrostatic interaction with the crystal lattice above. 

In order to understand the mechanisms of lattice excitation we reproduced the Mach 
cones with the molecular dynamics simulation code of Ref. 29[. It solved the equations of 
motion of 3000 microparticles (with the same mass as in the experiments) interacting with 
each other via a Yukawa potential. They were confined in a parabolic potential well, strongly 
in the vertical and weakly in the horizontal directions, forming a monolayer. Parameters 
of the simulation are listed in Table [H The lattice was excited by a 3D repulsive spot force 
field moving at 40 mm/s along the x-axis. The force field was similar to that of a charged 
particle. 

Both experimentally obtained (Fig. CLk-d) and simulated (Fig. [Tfe,f ) Mach cones have been 



visualised using the method of Refs. 
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, jjj • Several consecutive frames have been aligned at 
the cone vertex and averaged, producing high resolution velocity maps (Fig. [I]). This reduced 
the thermal noise and enhanced the motion correlated with the Mach cone. The horizontal 
speed is shown in Figs. [T^,c,e, it visualises compressional Mach cones (Figs. [T^,e). Note 
that there is no compressional cone in Fig. [lb, since the perturbing particle was subsonic. 
A double cone is visible in Fig. [T^,. The front compressional cone agreed well with the 
theoretical prediction (dashed line). The second cone with a smaller angle is a rarefactional 
cone [12 ]. 

Since the particles in a quasi-2D complex plasma are confined vertically in a finite har- 
monic potential, they can move in the vertical direction, e.g. out of the lattice plane. A 
fast charged particle moving under the main lattice layer exerts a 3D force and therefore 



perturbs the lattice vertica 
lations 27] and waves 



ly as well as horizontally. While large scale vertical lattice oscil- 



8[ have been observed before, here a long and narrow stripe of 
lattice is excited vertically. Figures [Tb,d,f show the vertical velocity maps of lattices in the 
wake of an out-of-plane perturbing particle. It is worth noting (to avoid misunderstanding) 
that a vertical velocity component mapping on the x, y plane is convenient way to reveal the 
vertical component of a Mach cone in a monolayer complex plasma. The maximal vertical 
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FIG. 1. Velocity maps visualising experiment I (a,b), experiment II (c,d), and simulation (e,f), 
see TabJT] for parameters. The maps show the absolute value of the particle horizontal velocity 
Vh (a,c,e) and the particle vertical velocity v z (b,d,f). The black dashed line in (a,b) represents 
the theoretical longitudinal Mach cone. The perturbation moved from left to right in all cases. A 
compressional Mach cone with multiple structure is clearly visible in cases, where the longitudinal 
Mach number was > 1 (a,e). The vertical wake structure (b,d,f) shows the vertical transverse waves 
generated at the cone vertex by the perturbation. The maps were obtained by aligning consecutive 
frames at the cone vertex and averaging them in order to reduce noise, 25 frames were averaged 
for experiment I (a,b), 13 frames - for experiment II (c,d), and 10 frames were used for simulation 
(e,f). 

particle displacements observed remained relatively small and, of course, not enough to re- 
construct precisely the vertical structure of a Mach cone (as done in the experiments with 



3D complex plasma clouds 



14jl). Both experiments and the simulation produce the same 



vertical excitation pattern, with a wavelength close to 3a and a width of the excited stripe 
of 3 — 5a where a is the average interparticle distance. This width remained nearly constant 
at different distances from the cone vertex indicating that the VTW speed is so low, that 
the vertical wave is damped by the neutral drag before it spreads noticeably. The theory- 
predicted VTW speed was a small percentage of the compressional dust-lattice wave speed. 
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Interestingly the vertical oscillations persisted further away from the cone vertex, than the 
compressional or shear waves, indicating that they experienced less damping (Fig. [Tb,d,f). 

A linear analytical theory has been used to describe the vertical lattice oscillations. Since 
the strength of vertical confinement of a stable monolayer is normally about 50 times higher 
than that of horizontal confinement, the amplitude of the vertical oscillations is much smaller 
than that of the horizontal oscillations, given that the energy is evenly distributed over the 



oscillation modes 



30j. The amplitude of the vertical motion is also much smaller than the 



lattice constant. The vertical displacement z A of the particles in a lattice obeys: 



d tt z A + ld t z A + [Q 2 ,z A - 2Q 2 C a A _ A >z A >] = 

A' 

= /(A-Vt) 



(1) 



where A = mai + n&2 {m and n are integers, ai^ are the primitive translation vectors 



3l|), the prime denotes the neighboring particles, 7 is the damping rate, fl v = 2nf v is the 

Z 2 e 2 

vertical confinement parameter, VL 2 C = -5- is the dust lattice frequency, a is the dispersion 



operator (with the Fourier spectrum «k (32 



33|). and / is the excitation force due to a 



repulsive projectile moving with a velocity V. The projectile is a point charge moving at a 

I ^ I 

height z p below the lattice and interacting with it via a force /k oc exp( 



A 



(kA)*). 




FIG. 2. 2D Fourier transform of the vertical velocity map obtained in experiment I, averaged over 
10 frames (a), and its theoretical prediction (b). 
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The solution of ([T]) is a convolution of the exciter Fourier spectrum / k , and the resonance 
propagator G k where the latter depends on 7: 



Zk = e 4kV 'G k / k , 
Gf k =[^-(kV + i 7 )kV]- 1 , 

nl = n 2 v - 2fi c 2 « k . 



(2) 
(3) 
(4) 



The Fourier-transformed spectral intensity of the vertical oscillations is, J k v kV k , v k = 
dtZ^. In our case of low damping, it shows a well-pronounced signature - a narrow stripe 
transverse to the direction of the projectile propagation in both experiment and theory 
(Figj2| in Figfj^a), scattered bright dots, not well aligned vertically, are most probably the 
traces of previous excitations with different directions of propagation). The position of the 
stripe corresponds to the resonance condition V p h = V, where V p ^, is the vertical wave 
phase velocity. Therefore the vertical structure of the Mach cone can be used to determine 
the resonant wave number. It was found to be ka = 2.02 ± 0.08 in the experiment, which 
agrees with the theoretical value of ka = 2.12. The spatial decrement of the resonant vertical 
oscillations ([2]) is L~ l ~ 7/V. Since the projectile velocity V is relatively large, the spatial 
decrement is low, and the wave pattern persists for a long distance ~ L explaining the 
seemingly lower vertical damping (in experiment I, L ~60 mm). 

In conclusion, we have reported the first direct observation of the vertical velocity com- 
ponent of a Mach cone in a quasi-2D complex plasma crystal excited by a particle moving 
underneath the main layer. The vertical structure or z-Mach cone was present in both cases 
of a super- or sub-sonic projectile. It is explained by excitation of the vertical transverse 
wave in the wake of the projectile. z-Mach cones can be used to determine the resonance 
wavelength of the vertical transverse wave. 
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